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ABSTRACT 

We derive expressions for the time integrated spectrum of Cosmic Rays (CRs) that are accelerated 
in a decelerating relativistic shock wave and escape ahead of the shock. It is assumed that at any given 
time the CRs have a power law form, carry a constant fraction of the energy E of the shocked plasma, 
and escape continuously at the maximal energy attainable. The spectrum of escaping particles is highly 
sensitive to the instantaneous spectral index due to the fact that the minimal energy, Emin ^ T^mpC^ 
where F is the shock Lorentz factor, changes with time. In particular, the escaping spectrum may be 
considerably harder than the canonical N{e) oc spectrum. For a shock expanding into a plasma 
of density n, a spectral break is expected at the maximal energy attainable at the transition to non 
relativistic velocities, e ~ 10^^(e_B/0.1)(n/lcm~^)^/^(£^/10^^erg)^/^eV where 63 is the fraction of the 
energy flux carried by the magnetic field. If ultra-high energy CRs are generated in decelerating 
relativistic blast waves arising from the explosion of stellar mass objects, their generation spectrum 
may therefore be different than the canonical N{e) oc e~^. 
Subject headings: cosmic rays 



1. INTRODUCTION 

Cosmic rays (CRs) are widely believed to be accel- 
erated in coUisionless shock waves in various sy stems 
(see e.g. IBlandford fc Eichleii IMtI: lAxfordl [1991 . In 
particular, supernova remnant (SNR) shocks propagat- 
ing subrelativistically are thought to accelerate pro- 
tons up to ~ 10^^ eV (see e.g. Gi nzburg fc Svrovatskiil 
IT969I) . while non or trans-relativistic internal shocks 
or ultra-relativistic external shocks of the jets of ac- 
tive galactic nuclei (AGNs, see e.g. [B crczinskv 2008) or 
gamma-ray bursts (GR Bs, iWaxman Ill995t iVietriill995i : 
IMilgrom fc "Usovl[T995l) may be the sources of Ultra High 
Energy CRs (UHECRs) up to the GZK energies of - 
few xlO^o eV. 

In cases where the shocked material expands consider- 
ably before releasing the CRs, most of the cosmic rays 
loose most of their energy by adiabatic losses before es- 
caping. Thus, the instantaneous spectrum of CRs at a 
given time and the integrated spectrum of the CRs es- 
caping from the system may be different. 

In this paper, we write down simple analytic expres- 
sions for the spectrum of escaping CRs from a relativistic 
decelerating shock wave assuming that at any given time 
(a) CRs escape at the maximal energy attainable at that 
time, (b) the CRs have a power law energy spectrum, 
and (c) the CRs carry a constant fraction of the shocked 
thermal plasma. We show that even under these simple, 
widely acceptable assumptions, the resulting spectrum 
may be non trivial and in particular significantly harder 
than the canonical N{e) oc spectrum. 

In section § [5] we formulate our basic assumptions and 
give general expressions for the resulting flux. In section 
§[3] we focus on energy conserving blast waves expanding 
in a uniform medium. The implications of the results 
are discussed in § |4l Expressions for the escaping CR 
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spectrum expected in other scenarios, including power 
law density profiles, radiative shocks, and a jet injection 
with constant luminosity, are derived in El 

2. BASIC MODEL AND ASSUMPTIONS 

Consider CRs that are accelerated by an expanding 
shock wave which is decelerating due to interaction with 
an external medium. 

2.1. Assumptions 

(i) The basic assumption that we make is that at any 
given time (or shock radius _R), particles escape at the 
maximal energy Emax attainable at that time. Specifi- 
cally, it is assumed that eiVesc(e), the number of CRs 
ejected within a logarithmic energy interval around e, 
is similar to the number of particles in the system, 
eN{e, -R), at the radius R for which the maximal energy 
£max is equal to e. 



£7Vesc(£) -^eiV(e,i?|e„..=s). 



(1) 



Here and everywhere else in this letter, all quantities 
are measured in the observer frame. This assumption 
is expected to be true whenever the acceleration is lim- 
ited by the finit e size of the system (for a recent dis- 
cussion see e.g. iCaprioh et"alTl2009D . Note that for 
a relativistic shock with Lorenz factor F, only parti- 
cles that move in the shock propagation direction to 
within an angular deviation of 1/F, can outrun the shock 
and escape. This does not introduce a further correc- 
tion however, due to the fact that all particles in the 
shocked region, including the thermal and accelerated 
particles, are beamed in the observer frame to an an- 
gular separation of 1/F (for examples of the expected 
ang ular di stribution of CRs in relativis tic shocks see 
e.g. iKirk fc Schneider 1987; Bcdnarz & Ostrowskil 11991 
iKirk et al.ri2000; .Keshet fc Waxman.2005il . 

For concreteness we further make the following as- 
sumptions. 

(ii) At any given radius R, the energy spectrum of CRs 
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is a power law, N{e) oc e 



-2-x 



for £,„in < £ < 



with 



a; > and with total energy Ecr- The spectrum can be 
expressed as 



e^N{e) = xEcR 



where we assumed that (emax/Smin) ^ ^ 1- 



(2) 



(iii) The minimal, maximal and total cosmic ray energies 
are power law functions of the radius, 

e„i„ cx i?-"""", £^ax«i?-"— , EcR^R-"^. (3) 

2.2. Resulting spectra 

Under the above assumptions, the spectrum of escaped 
particles is given by 

e^N,,,ie) oc e""-" (4) 



with 



= a; — (a-a 



(5) 



For the case a^,i„ = 0, this reduces to equation (28) of 
lOhira et al.l (|200l . 

Equation ^ is valid for a; > 0. We note that for the 
limiting value, a; = 0, a logarithmic correction is intro- 
duced to the spectrum of escaped CRs, due to the loga- 
rithmic dependence of the total CR energy on £„ 



-min,maxi 



EcR = log(£inax/emin)e^(e)- For a; < 0, the energy 
carried by the CRs is dominated by the particles with 
largest energies, e = Smax- The resulting escaped spec- 
trum satisfies a;osc — o:e /ctmax and is not sensitive to the 
form of the instantaneous spectrum. 

3. ENERGY CONSERVING BLAST WAVES 

We next apply the above general result to simple cases, 
focusing on a blast wave of fixed energy which starts off 
ultra-relativistic, and then decelerates to non relativis- 
tic velocity. The non relativistic and relativistic stages 
are analyzed in § 13.11 and § 13.21 respectively. The spec- 
trum resulting from the combination of the two stages is 
addressed in § SI 

3.1. Non-relativistic expansion 

Perhaps the simplest case to consider is the case of 
constant CR minimal energy £min and total energy Eqr 
(aE,amin = 0). This is expected in the Sedov- Taylor 
phase of non relativistic blast waves. The energy in CRs 
is dominated by the relativistic CRs, e > mpC^, and is 
commonly assumed to be a constant fraction / of the 
total energy E in the system. Hence Emin ~ m^c^ and 
-E-CR = fE, both independent of the shock radius. Eq. 
([U implies that the escaped spectrum is the same as the 
instantaneous spectrum. 



I.e. iVesc(£) « £-2-.^ (g) 

This result can be obtained directly from Eq. (IT|) by not- 
ing that the entire spectrum, up to the maximal energy 
emax(^), is independent of radius. 

3.2. Ultra-relativistic expansion 

We next consider the case of ultra-relativistic expan- 
sion. Below and in § |X] we make the following assump- 
tions in addition to assumptions (i)-(iii) above: 



(iv) The minimal CR energy is the energy of the thermal 
particles, 

emin ~ £th ~ r^mpC^; (7) 

(v) The CR pressure is a radius independent fraction fcR 
of the momentum flux in the shock frame, pcR oc T^p, 
implying 

Ecr cx R^'V^p; (8) 

(vi) The maximal CR energy is the maximal energy of 
CRs that are confined by a fluid rest frame magnetic field 
Srost (equivalent to Diffusive Shock Acceleration in the 
Bohm limit), 

Emax - eEobsSR - eB,,,tR « eiTp^/^)R, (9) 

where -Bobs ^ FSicst is the electric field in the ob- 
server frame corresponding to a magnetic field -Brest ~ 
{SnT'^ pc^eBY^"^ in the rest frame of the shocked plasma, 
assumed to carry a constant fraction es of the momen- 
tum flux, and SR ^ R/T is the maximal distance that a 
cosmic ray can propagate along the electric field in the 
observer frame. 

3.2.1. Ultra-relativistic impulsive expansion 

We next consider a Blandford- McKee (energy conserv- 
ing; Blandford & McKec 1976) shock expanding into a 
uniform medium. By assumption = 0. Using equa- 
tions d!])-®, c^max — 1/2 and ctmin — 3 are obtained. 
Using ([S]) we find a;osc = —5a:, or 



iVosc(£) oc £ 



bx~2 



(10) 



4. DISCUSSION 

The spectral index of escaped CRs, given by Eq. (fTO| . 
may be surprisingly hard if the instantaneous spectrum 
(equ.[2]) is softer than a fiat spectrum, i.e. a: > 0. The 
basic reason for this is the fact that the minimal CR 
energy ~ T^rUpC^ changes with radius much faster 
than the maximal energy, implying that at later times 
(corresponding to lower escaped energies) there is much 
less CR energy at the escaping, high end of the spectrum. 

For example, consider the commonly a ssumed Diffusive 
Shock Accele r ation (DSA) mech a nism (iKrvmskii 19771: 
Axford et all [19771 [Belli [l978l : IBlandford fc Ostriked 



19781), which for ultra-relativistic shocks with isotropic, 
small-angle scattering, leads to an instantaneous 
spectrum N oc £~'^°^^_ (iKeshet fc WaxmanI 120051 : 
iBednarz fc Ostrowskil[T99l iKirk et al.ll2000f ). Using Eq. 
(fTO|) , this implies a very hard escaping spectrum, Ncsc oc 
£-8/9^ We emphasize that DSA has not been shown to 
work based on first principles and even if it does, the 
resulting spectrum is se nsitive to t he scattering mecha- 
nism (e.g. iKeshe t fc WaxmanI 120051 : [keshet 20061 which 
is poorly understood. In fact, the instantaneous spec- 
trum doesjio^jiecessaxiIyneed_t^^ form 
(e.g. iCaprioH et alll2009l: lOhira et al.ll2009D or to h ave a 
constant spectral index(e.g. lEllison fc Doubl^ l20041. 
Note t hat an instantane ous spectrum th at is fi at N oc 
(e.g. iKatz et al.l [2007l) or harder fe.g. lKeshetll200ll : 
IStecker et al.ll2007r and others) can be obtained in spe- 
cific particle acceleration models. These models, how- 
ever, are different from the possible hard spectrum of 
escaping particles discussed above. In fact, such hard 



CR spectrum 
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instantaneous spectra lead to a flat escaped spectrum 
Ncscis) fx [see discussion following eq. ([5])]. The 
important point is that the escaping spectrum is highly 
uncertain, very sensitive to the acceleration mechanism, 
and may be considerably harder than Niseis) oc 

Once the shock decelerates to non relativistic speeds, 
the integrated escaping spectrum changes its form. In 
the subsequent Sedov- Taylor phase, the minimal CR en- 
ergy does not change any more, Smin ~ mpC^, and the 
escaping spectrum will be similar to the instantaneous 
spectrum [see Eq. ([5])]. Note that in the extreme case 
in which CRs carry most of the energy, and assuming 
they are accelerated by non- lin ear DSA (for a recent re- 
view see lMalkov fc Dnirvll200H ). the instantaneous spec- 
trum may not be a power law and correspondingly the 
escaping spectrum i s different th an that given here (e.g. 
ICaprioli et all 120091: lOhira et al. 'WOQl. In any case, a 
significant break is expected in the escaping spectrum at 
an energy equal to the maximal energy achievable at 
the transition from relativistic to non-relativistic shock 
velocities. Using Eq. (|9]) for r/3 ~ 1, the break is ex- 
pected at 

^ lO'^eB^^ml/'EU'eY, (II) 



where we have taken, as an example, a gamma-ray burst 
source with E = 10''^"'^ £'51 erg, p = n rup is the ambient 
density with n = uq cm~^, and es = O.Ie^.-i. The spec- 
tral index above this energy will be —2 — Xcsc — ~2 -I- 5x, 
as given by Eq. (ITU)) . We note that the energy range 
where Xcsc — ~5x is limited between and max(£) ~ 
emax(ro) where To is the initial Lorentz factor. The max- 
imal energy scales with T as £max F^/^ resulting in a 

1/3 

narrow energy range e, < e < Eg' While this range 
depends on model parameters, it does point to a possible 
mechanism for achieving a hardening of the spectrum at 
the high end of the CR spectrum. 

We conclude that the spectrum of UHECRs, may be 
different than e^^, and possibly considerably harder, e.g. 
if they originate from a relativistic decelerating blast- 
wave resulting from the explosion of a stellar mass object. 
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TABLE 1 

Spectral indexes of escaping CRs for different relativistic scenarios 



Scenario 


Conserved quantity 


dp 


Xcsc [p oc _R"] 


Xcsc [poc R 2] 


Impulsive, Energy conserving 


E ~ r'^pR^ 


(3-ap)/2 


—5a; 


— X 


Impulsive, Radiative 


MT ~ FpR^ 


3 - ap 


-3/2 - 2x 


-1/2 - X 


Continuous injection L=const. 




(2-ap)/4 


-4+3x 


2-x 



The value of the spectral index Xesc is given by Eqs. IIASI l and (IA6I I for the various scenarios. N{e) oc e ^ 

APPENDIX 
OTHER PARTICULAR CASES 
Consider an ultra-relativistic shock, which expands into a medium with a density profile 

p oc R-"" 

with a bulk Lorentz factor 



r oc 



Using Eqs. dZ])-® and (|XT| we find: 
Substituting this in equation ([5]) we find 



3 — 2ar — ctn 



1, and CUE = 2ar + — 3. 



(2-ap)+2ar 



(2-a,)-2ar-^- 

For a uniform density distribution, ap = 0, and for a wind-like density distribution, = 2, this reduces to 

3 — 2ar 1 + ar 



and 



respectively. The escaping spectrum is then 



ar — 1 i — ar 

_ 1 - 2ar 

'^csc \ Z~ ^ 

ar + 1 

A^e.c(e) oce-2-Xe.c^ 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 

(A6) 
(A7) 



which is flatter then —2 for Xesc < 0. The expected spectral indices x^sc, for various physical scenarios, are provided 
in table [T] 



